Geometric Convergence of the Metropolis Hastings Simulation Algorithm Lars Holden Norwegian Computing Center and University of Oslo
نویسنده
چکیده
Necessary and su cient conditions for geometric con vergence in the relative supremum norm of the Metropolis Hastings simulation algorithm with a general generating function are estab lished An explicit expression for the convergence rate is given Introduction This paper discusses the convergence rate for the Metropolis Hastings simulation algorithm proposed in Hastings The Metropolis Hastings simulation algorithm is used for sampling from a distribution f x There is currently a lot of interest in MCMC both theoretically and in a large number of applications see Geyer The challenge in Metropolis Hastings is to nd a good generating function The explicit formula for the convergence rate given in this paper may be used to compare di erent generating functions Meyn Tweedie prove that the Doeblin condition is equivalent to uniform ergodic i e uniform convergence in total Date April Mathematics Subject Classi cation J
منابع مشابه
Approximating Bayes Estimates by Means of the Tierney Kadane, Importance Sampling and Metropolis-Hastings within Gibbs Methods in the Poisson-Exponential Distribution: A Comparative Study
Here, we work on the problem of point estimation of the parameters of the Poisson-exponential distribution through the Bayesian and maximum likelihood methods based on complete samples. The point Bayes estimates under the symmetric squared error loss (SEL) function are approximated using three methods, namely the Tierney Kadane approximation method, the importance sampling method and the Metrop...
متن کاملImproving convergence of the Hastings - MetropolisAlgorithm with a learning proposal
The Hastings-Metropolis algorithm is a general MCMC method for sampling from a density known up to a constant. Geometric convergence of this algorithm has been proved under conditions relative to the instrumental distribution (or proposal). We present an inhomogeneous Hastings-Metropolis algorithm for which the proposal density approximates the target density, as the number of iterations increa...
متن کاملAdaptive Independent Metropolis – Hastings
We propose an adaptive independent Metropolis–Hastings algorithm with the ability to learn from all previous proposals in the chain except the current location. It is an extension of the independent Metropolis–Hastings algorithm. Convergence is proved provided a strong Doeblin condition is satisfied, which essentially requires that all the proposal functions have uniformly heavier tails than th...
متن کاملOperator splitting methods for degenerate convection–diffusion equations II: numerical examples with emphasis on reservoir simulation and sedimentation
Helge Holden a, Kenneth Hvistendahl Karlsen b and Knut-Andreas Lie c a Department of Mathematical Sciences, Norwegian University of Science and Technology, N-7491 Trondheim, Norway E-mail: [email protected] b Department of Mathematics, University of Bergen, Johs. Brunsgt. 12, N-5008 Bergen, Norway E-mail: [email protected] c Department of Informatics, University of Oslo, P.O. Box 1080 Blinde...
متن کاملComputing the Matrix Geometric Mean of Two HPD Matrices: A Stable Iterative Method
A new iteration scheme for computing the sign of a matrix which has no pure imaginary eigenvalues is presented. Then, by applying a well-known identity in matrix functions theory, an algorithm for computing the geometric mean of two Hermitian positive definite matrices is constructed. Moreover, another efficient algorithm for this purpose is derived free from the computation of principal matrix...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 1998